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We calculate the spin relaxation rates in a parabolic InSb quantum dots due to the spin in- 
teraction with acoustical phonons. We considered the deformation potential mechanism as the 
dominant electron-phonon coupling in the Pavlov-Firsov spin-phonon Hamiltonian. By studying 
suitable choices of magnetic field and lateral dot size, we determine regions where the spin relax- 
ation rates can be practically suppressed. We analyze the behavior of the spin relaxation rates as a 
function of an external magnetic field and mean quantum dot radius. Effects of the spin admixture 
due to Dresselhaus contribution to spin-orbit interaction are also discussed. 

PACS numbers: PACS: 73.20.Dx, 63.20.Kr, 71.38.+i 



I. INTRODUCTION 



The ability to manipulate and control processes that 
involve transitions between spin states is, at the moment, 
of extreme importance due to the recent applications 
in polarized spin electronics and quantum computation. 
Spin dephasing is the most critical aspect that should 
be considered in the elaboration of proposals of quan- 
tum computation based in single spin states as qubits in 
quantum dots (QDs) Currently, QDs of diverse ge- 
ometries are considered as promising candidates for im- 
plementation of quantum computation devices because 
the electronic, magnetic and optical properties can be 
controlled through the modern grown and nanofabrica- 
tion techniques. Due to the long electron spin dephas- 
ing time experimentally reported 0, S 0] ) the spin of an 
electron localized in a QD has been suggested to realize a 
quantum bit. While for bulk and for 2D systems the spin 
relaxation processes have been studied in some detail, the 
problem for QD's still require deeper and further discus- 
sions. Several processes that can induce spin relaxation 
in semiconductors have been identified and were stud- 
ied. At the moment remains in discussion which, between 
these processes, is dominant in zero-dimensional systems. 
Some experimental results have shown good agreement 
with the theoretical predictions for 2D systems |5| but, in 
general, the identification of the processes through direct 
comparison with the experimental results may become a 
formidable task. This problem is more critical for QDs, 
since few experimental results exist and the theoretical 
discussion of the spin relaxation mechanisms is still an 



open subject. Extensive theoretical works in QD sys- 
tems have studied the main phonon mediated spin-flip 
mechanisms, including admixture processes due to spin- 
orbit coupling and phonon coupling due to interface 
motion (ripple mechanism) Q . Spin relaxation rates de- 
pend strongly on the dot size, magnetic field strength, 
and temperature, as reported by several authors [gl, [1]. 
It was shown that the quantum confinement produces, in 
general, a strong reduction of the QD relaxation rates. 

In this work, we calculate the spin-flip transition rates, 
considering the phonon modulation by the spin-orbit 
interaction. For this purpose we will use the spin- 
phonon interaction Hamiltonian proposed by Pavlov and 
Firsov [13] ■ In this model, the Hamiltonian describing 
the transitions with spin reversal, due to the scattering 
of electrons by phonons, can be written in a general form 
as 



Hph = Vpi, + -f[ax VVph] • (p + e/cA), 



(1) 
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where V^h is the phonon operator, a is the Pauli spin 
operator, 7 is related with the strength of the electron- 
phonon interaction, p is the linear momentum operator 
and A is the vectorial potential related with the external 
magnetic field B. This model has the advantage of being 
easily adapted to the study of other interaction mecha- 
nisms with phonons. 



II. THEORY 

Based on the effective mass theory applied to the prob- 
lem of the interaction of an electron with lattice vibra- 
tions, including the spin-orbit interaction and in presence 
of an external magnetic field, Pavlov and Firsov [3, IToj 
have obtained the spin-phonon Hamiltonian that de- 
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scribes the transitions with spin reversal of the conduc- 
tion band electrons due to scattering with longitudinal 



J 



lattice vibrations as 
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where, &q(&q) are annihilation (creation) phonon opera- 
tors, the magnetic vector potential A is obtained in the 
symmetric gauge considering an external magnetic field 
B oriented along the z axis, = x ± iy, where x, y 
are unitary vectors along the x and y axis. Cq is a unit 
vector in the direction of the phonon polarization, q is 
the phonon wave vector, p is the momentum operator, v 
is the average sound velocity, pM is the mass density, V 
is the system volume and d{q) is a deformation potential 
electron-phonon coupling constant [To| 



can be written as 
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where i?g=0.17eV is de band gap for InSb, TO*=0.013mo, 
Ea—7.0 eV is the deformation potential, A=0.81 eV is 
the spin-orbit splitting of the valence band, and f{x) — 
(l + 2a;)/[(l + |x)(l+x)]. 

According to Eq. ([2]), we will restrict ourselves to the 
case of bulk phonon modes. The longitudinal acous- 
tic phonons are described by plane-waves of the form 
UoBq exp*'! '', where Uq is the normalization constant. We 
consider the phonon energy dispersion in the Debye ap- 
proximation: Eq — hvq, for InSb we use v = 3.4 x 10'^ 
m/s. It is important to point out that the theoretical 
treatment of the spin-phonon interaction developed by 
Pavlov and Firsov 0, can be easily adapted to sev- 
eral electron-phonon interaction processes. An appropri- 
ate choice of the electron-phonon coupling d{q) allows 
us to adapt the matrix element ([2]) for piezoacoustic or 
LO-optical spin scattering. 

It has been assumed that the confinement along the z 
axis is much stronger than the lateral confinement. Thus, 
the lateral motion is decoupled from the one along z and 
the envelope functions separate ipir) — f{x,y)(l){z). The 
^-dependent part of ?/'(r) is an eigenfunction of a sym- 
metric quantum well of width L. In lens-shaped quasi- 
two dimensional self assembled QDs, the bound states of 
both electrons and valence-band holes can be understood 
by assuming a lateral spatial confinement modeled by a 
parabolic potential with rotational symmetry in the x — y 
plane jTl|, V{p) = imwoP^, where huiQ is the character- 
istic confinement energy, and p is the radial coordinate. 
By using the one-band effective mass approximation and 
considering an external magnetic field B applied normal 
to the plane of the QD, the electron lateral wave function 
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where Cn,M = \/n\/['K{n -f lJ^'^' is the Laguerre 

polynomial, n (M) is the principal (aziniuthal) quantum 
number, and xi'^) is the spin wave function for the spin 
variable a. The corresponding eigenenergies are 

En M a = {'2n+ \M\ + l)hn + {M/2)hwc + {cr/2)gpBB, 

(5) 

where i7 ~ (cjq -t- uj^/^y^^, /is is the Bohr magneton, 
a = {h/mfiy/^ is the effective length, a — ±1 for spin 
up and down respectively, and tOc — eB/m. 

In our model, we also consider the effects of the Dres- 
selhaus contribution that provides additional admixture 
between spin states. For 2D systems, the linear Dressel- 
haus Hamiltonian can be written as 



Bd = [CfxPx - O-yPy) : 



(6) 



where pi — ~ih\J i ~\- {e/c)Ai with i — x,y and (3 is the 
Dresselhaus coupling parameter for this confinement. If 
the confinement potential in the z-direction is considered 
highly symmetrical, then VV^ ^ and the Rashba con- 
tribution can be safely ignored. 

The spin relaxation rates (W) between the electronic 
states: (rt, M, f (|)) (n',Af',J, (|)), with emission 
of one acoustic phonon, arc calculated from the Fermi 
golden rule. 

In the Hamiltonian ^ , we only consider the deforma- 
tion potential electron-phonon coupling, this is due to the 
large g-factor in narrow gap InSb (\g\ ~ 51), the domi- 
nant electron-phonon coupling for spin relaxation is the 
deformation potential mechanism [12| . The piezoelectric 
coupling governs the spin relaxation processes in wide or 
intermediate gap semiconductors. In the transition ma- 
trix elements calculation, we not only consider the linear 
term iq • r in the expansion of exp(iq • r) '6], but also 
the integral representation of Bessel function is used in 
the evaluation of electron-phonon overlap integrals. The 
linear approximation of exp(iq • r) may be valid for spin 
inversion transitions in the spin polarized ground-states 
of GaAs based QDs where, due to the small value of the 
electron factor, only long wavelength phonons are in- 
volved. 
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III. RESULTS AND DISCUSSION 

The calculations were performed for a parabolic InSb 
QD at T ^ K. The material parameters for the InSb 
system are listed in Ref. [11]. We only have considered 
electron transitions between ground state electron Zee- 
man levels (0,0, T) ^ (0,0, J.) and (0,1, j) ^ (0,1, T). 
The temperature dependence for one-phonon emission 
rate is determined from W = Wo{nB + 1), where ub 
is the Bose-Einstein distribution function and Wq is the 
rate at T = K. In the temperature regime T <10 K, 
we obtain + 1 « 1 and W ~ Wq. For temperatures 
larger than few Kelvin degrees (T > 10 K), two-phonon 
processes become the dominant spin relaxation mecha- 
nism Q • These types of processes have not been consid- 
ered in the present calculation. 
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FIG. 1: (Color online) Spin relaxation rates, W, for a 
parabolic InSb QD considering the deformation potential cou- 
pling mechanism as a function of the magnetic field B. We 
consider the transitions: (0, 0, |) (0, 0, i) and (0,1,1) 
(0, 1, t) and several lateral dot radius ro = 50, 75, 100, 125, 
150, 175, and 200 A(same ro ordering for both transitions) 

In the Fig. [T] we show the spin relaxation rates due to 
deformation potential electron-phonon mechanism, as a 
function of the external magnetic field B and considering 
so me typi cal values for the effective lateral QD size, rg — 
a/ h/mujQ. 

Some interesting facts about these results should be 
pointed out: 

i) The rates show a strong dependence with the mag- 
netic field. This fact can be explained from the depen- 
dence of the rates with the transition energy AE. In gen- 
eral, we obtain that W ^ [g*fiBB]^ — (AE)", n being an 
integer number that depends on the electron-phonon cou- 
pling process and g* the effective (7-factor. As can be seen 
in Fig. [1] when the magnetic field increases, the rates 
also increase until reaching a maximum near B ~ 0.5 T. 
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FIG. 2: (Color online) Contour plot of the spin relaxation 
rates as a function of magnetic field B and lateral size ro for 
transitions: a) (0, 0, T) ^ (0, 0, [) and b)(0, 1, ],) ^ (0, 1, T). 



The position of this main maximum is defined from the 
transition energy conservation: Enia' — En'i'a-' = hvq. 

ii) The oscillatory behavior of the rates, observed for 
B > 0.7T are mainly produced by complex interplay be- 
tween the spin admixture and electron-phonon overlap 
integrals. The phonon modulated Rashba interaction, 
given in Eqs. (H]) and ^ produce level admixture ac- 
cording to the selection rules: (n, M, f) — > {n, M + 1, |), 
(n,M,i) ^ (n,M-l,t) and (n,M,T (i)) ^ (n, M, i 
(t)). The last condition, which hybridize states with 
different spin orientations at the same level M, is the 
main responsible for the rates oscillations at _B > 0.7T. 
For AM — transitions, the Dresselhaus SO interaction, 
given in ©, is not able to produce spin admixture. 

As it is shown in Fig. [1] the (7*-factor effects are par- 
ticulary important for the ground-state Zeeman transi- 
tion. For small magnetic fields, g* — > (?buik and we may 
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neglect the spin admixture effects. Therefore, the spin 
relaxation shows no oscillations and becomes almost in- 
dependent of tq. This small QD size dependence is in 
agreement with the experimental observations of Gupta 
and Kikkawa fll|. 

iii) The rates dependence with the lateral QD size rg, 
are related to the interplay effects between the spatial 
and magnetic confinements. This competing effects are 
contained in the electron-phonon overlap integral, / cx 
/ /n',i',a'(p)exp(«q ■ r)fn,i,a{p)dr. For large fields, the 
magnetic confinement causes a gradual decrease in the 
overlap integral as the tq increases. For small magnetic 
fields, the spatial confinement is dominant. Thus, when 
ro diminishes the wave functions become more localized 
and the overlap integral should increase. This effects 
explain the behavior of the spin transition (0,1, |) — > 
(0, 1, t) showed in Fig. [1] (red lines). The Zeeman ground- 
state rates (black lines) are strongly dependent on AE 
and, for small B, the rates are weakly dependent on /. 

iv) The same rates calculated for GaAs (not showed 
here), are in general, one order of magnitude smaller than 
InSb rates. As we expected, the relaxation via piezoelec- 
tric coupling is more efficient than via the deformation 
potential phonon processes. 

In Fig. [2] a) we have plotted the spin relaxation rates 
for the ground-state Zeeman transition as a function of 
To and B. We clearly identify a region of strong spin 
coherence, defined by S > 1 T and tq > 100 A. In this 
regime, the relaxation times r arc in the ns order and this 
is an important feature for spin qubit engineering. In the 
B < O.IT regime, the relaxation times are approximately 
of few /is. This spin frozen region is not robust against 
the temperature and will disappear whenever the thermal 
energy is larger than the spin transition energy. The plot 
in Fig. Ob) shows the spin rates for (0, 1, |) — > (0, 1, t) 
transition. As in the previous case, the strong coherence 



regime is defined approximately by B > 1 T and tq > 100 
A. However, the relaxation times, in this region, are two 
orders of magnitude faster than for (0,0, t) (0,0, |) 
transition. 

In conclusion, in this work we presented the calculation 
of spin relaxation rates via electron-phonon interaction. 
The model used in our work, allow us to estimate the 
magnitude of the scattering rates considering the most 
important contributions to the spin relaxation through 
acoustic deformation potential interaction. The rates ex- 
hibit a strong dependence with the magnetic field and 
with the lateral QD size. We identify regimes of strong 
spin coherence, for magnetic field, i? > IT, and for lat- 
eral sizes, ro > 100 A, we obtain relaxation times ap- 
proximately of 10~^ ns to 1.0 ns. The Zeeman ground 
state transition (0,0, t) — > (0,0, |), exhibits longer spin 
coherence times than the high energy transitions, which 
is important for spin qubit applications. The magnitudes 
of the calculated rates are compatible with those ob- 
tained experimentally. For a more rigorous analysis of 
the spin relaxation process, many other effects should be 
considered, like the anisotropy of the electron g-factor 
and the effects of the spatial confinement on the acous- 
tic phonons. These two issues, are still under intense 
discussion and should be considered for a more precise 
determination of the spin relaxation times. 
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